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Cutoff-effects in the spectrum of dynamical Wilson fermions 

M. Delia Morte a , R. Hoffmann", F. Knechtli a , and U. WolfP 

a Institut fiir Physik, Humboldt Universitat, Newtonstr. 15, 12489 Berlin, Germany 

We investigate the low-lying eigenvalues of the improved Wilson-Dirac operator in the Schrodinger functional 
with two dynamical quark flavors. At a lattice spacing of approximately 0.1 fm we find more very small eigenvalues 
than in the quenched case. These cause problems with HMC-type algorithms and in the evaluation of fermionic 
correlation functions. Through a simulation at a finer lattice spacing we are able to establish their nature as 
cutoff-effects. 



1. Motivation 

Recently more and more evidence has been 
accumulated that for dynamical improved Wil- 
son fermions at a lattice spacing of a ~ O.lfm 
the cutoff-effects are much larger than expected 
from quenched experience (for a summary of large 
scaling violations in the two-flavor-theory see 
ref. PP). As an extreme example, for three fla- 
vors the existence of a phase transition in the (3- 
K-plane has been numerically conjectured and is 
interpreted as a lattice artifact 

In addition several collaborations have re- 
ported algorithmic difficulties with (improved) 
Wilson fermions on relatively coarse lattices (see 
e.g. jj]), which seem to be related to small eigen- 
values of the Dirac operator. Through the inver- 
sion of the Dirac operator these small eigenvalues 
can result in large driving force during the molec- 
ular dynamics evolution of Hybrid Monte Carlo 
(HMC) algorithms. In turn, the large forces are 
likely to trigger instabilities of the numeric inte- 
grator employed, which produces large Hamilto- 
nian violations and can result in long periods of 
rejection 0]. In this way small eigenvalues affect 
algorithmic performance. 

The occurrence of very small eigenvalues gen- 
erates not only algorithmic problems. On the 
configurations in question the quark propaga- 
tor becomes large and one observes "spikes" in 
fermionic correlation functions. This affects not 
only their mean value but also their autocorrela- 
tion function, making the statistical analysis dif- 
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ficult. This is discussed in more detail in ref. @|. 

For dynamical simulations this is an unex- 
pected problem since naively it is assumed that 
the determinant should suppress small eigenval- 
ues compared to the quenched situation. In our 
setup we have two infrared cutoffs (finite quark 
mass and Schrodinger functional boundary con- 
ditions) that should prevent the Dirac operator 
from developing very small eigenvalues. Never- 
theless wc observe them in simulations with a 
lattice spacing of about O.lfm. Here we show 
that these small eigenvalues disappear if one goes 
to smaller lattice spacings and can thus be inter- 
preted as a lattice artifact. 

2. Setup and error analysis 

We simulate the Schrodinger functional (SF) 
with two dynamical flavors of non-perturbatively 
improved Wilson fermions. The algorithms used 
are HMC with two pseudo-fermion fields [3] and 
PHMC UJ. In the following the term 'eigen- 
value' always refers to the eigenvalue (in lattice 
units) of the square of the Hermitian even-odd- 
preconditioned Wilson-Dirac operator in the 
normalization of 0. For PHMC the parameters 
of the polynomial are chosen such that more con- 
figurations with small eigenvalues are produced 
compared to the QCD Boltzmann weight. After 
reweighting this gives a very good estimate of 
the path-integral weight of such configurations. 
Another benefit of using PHMC is that the poly- 
nomial provides a regularized inversion, thus also 
addressing the algorithmic problems due to large 
forces mentioned above. With this algorithm the 
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correct ensemble average of an estimator O is 
given by 
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where W is the reweighting factor [S] and the sub- 
script P indicates an average over the PHMC- 
generated ensemble. HMC-type algorithms are 
expensive and generate strongly autocorrelated 
data, which makes a careful data analysis indis- 
pensible. We use an explicit integration of the 
autocorrelation function as described in ref. [Hj- 

To obtain a histogram for a quantity / = (</>) we 
analyze P n — (x n (<?!>)}, where \n is the characteris- 
tic function of the n-th bin. This also provides us 
with an error on the population of the bin, thus 
giving us the ability to assess differences between 
histograms. If the PHMC algorithm is used, the 
bin population is given by 
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3. Comparison to the quenched case 



In order to test the naive expectation that the 
fcrmionic determinant suppresses small eigenval- 
ues we compare quenched and dynamical ensem- 
bles of 8 3 xl8 lattices at roughly matched physical 
parameters. Using the quenched data from ref. |5] 
and the dynamical data from refs. JU| an d d] we 
cohose the parameters such that the lattice spac- 
ing and the (large volume) pseudo-scalar mass 
are matched. 

The distribution of the smallest eigenvalue A m i n 
is shown in Figure 2] The quite small error bars 
we obtain from eq.[21at the lower end of the spec- 
trum are due to the enhanced occurrence of small 
eigenvalues by PHMC. While (A m i n ) is increased 
from 1.44(1) -lCT 4 to 1.72(5) -lO" 4 with two dy- 
namical flavors we see that in the infrared tail 
the dynamical data show more events. More pre- 
cisely, the probability of finding a smallest eigen- 
value below 4-10 -5 increases from 0.81(16)% to 
1.88(26)%. To show that this long tail towards 
zero is a cutoff-effect we now compare dynamical 
data from different lattice spacings at matched 
parameters. 
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Figure 1. Distributions of A m i n from matched 
quenched ((3 = 6.0) and two-flavor dynamical sim- 
ulations (/3 = 5.2). . 
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Figure 2. Distributions of the smallest eigenvalue 
from two-flavor simulations at two different lat- 
tice spacings. 

4. Finer lattices 

Through measurements of the SF coupling g 2 
|12I13) we found that increasing f3 from 5.2 to 5.5 
changes the lattice spacing by roughly a factor of 
2/3. Anticipating that the algorithmic problems 
due to very small eigenvalues would no longer be 
present at this finer lattice spacings we used HMC 
with two pseudo-fermions to generate an ensem- 
ble of 12 3 x 27 lattices. Ignoring small changes 
in the renormalization factor we compare this to 
an 8 3 x 18 PHMC ensemble at = 5.2, which is 
matched using LmpcAd the box size times the 
bare PCAC mass. To plot both lattice spacings 
simultaneously we divide by (A m ; n ) in Figure |21 

Going from j3 = 5.2 to 5.5 reduces the vari- 
ance (normalized by the mean value) of the small- 
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Figure 3. Comparison of the smallest eigenvalue 
distribution at a lattice spacing of approximately 
0.07 fm. 

est eigenvalue from 0.178(10) to 0.127(19). The 
long tail at the infrared has disappeared at the 
finer lattice spacing and we thus interpret it as a 
cutoff-effect. 

In a final step we compare the dynamical data 
at [3 = 5.5 to another quenched run at approxi- 
mately the same lattice spacing (/? = 6.26), vol- 
ume and bare quark mass. 

In Figure[3]both the quenched and the dynam- 
ical data show very similar behavior at the in- 
frared end of the spectrum. A comparison to Fig- 
ure n shows that the excess of very small eigen- 
values we found at the coarser lattice spacing has 
disappeared entirely. As at a ~ 0.1 fm the aver- 
age smallest eigenvalue is still shifted upwards by 
the fcrmionic determinant. 

5. Conclusions 

The problems in simulating dynamical Wilson 
fermions at a lattice spacing of approximately 
0.1 fm are due to the occurrence of very small 
eigenvalues in the spectrum of the Wilson-Dirac 
operator. We use PHMC to better sample this 
part of the spectrum in a comparison of two- 
flavor and quenched simulations at matched phys- 
ical parameters. 

As expected we find the average smallest eigen- 
value to be larger in the dynamical case. How- 
ever, due to its increased variance the dynamical 
data shows more very small eigenvalues despite 
the finite quark mass and the cutoff provided by 
the Dirichlet boundary conditions. 

At a lattice spacing of approximately 0.07 fm 
both the dynamical and a matched quenched run 



show a spectrum that is well separated from zero. 
As before (A m i n ) is larger for two flavors. 

We conclude that at (3 = 5.2, corresponding to 
a lattice spacing of 0.1 fm, the spectrum of the 
Wilson-Dirac operator is strongly distorted. Our 
simulations show that this is a cutoff-effect that 
disappears rapidly with increasing (3. 

We do not expect that these findings are spe- 
cific to the Schrodinger functional. Without the 
additional infrared cutoff due to the boundary 
conditions one should see these problems already 
at larger quark masses. 
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